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Abstract 

We propose and study a new first order version of the ghost-free massive gravity. Instead of 
metrics or tetrads, it uses a connection together with Plebanski's chiral 2-forms as fundamental 
variables, rendering the phase space structure similar to that of SU (2) gauge theories. The chiral 
description simplifies computations of the constraint algebra, and allows us to perform the complete 
canonical analysis of the system. In particular, we explicitly compute the secondary constraint and 
carry out the stabilization procedure, thus proving that in general the theory propagates 7 degrees 
of freedom, consistently with previous claims. Finally, we point out that the description in terms 
of 2-forms opens the door to an infinite class of ghost-free massive bi-gravity actions. 
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1 Introduction 



A deceptively simple way to modify general relativity in the infrared is to give the graviton a small 
mass. Massive gravity is a subject with a long history, and a most recent wave of interest comes 
from a possibility that it may provide a technically natural solution to the cosmological constant 
problem, see e.g. [T] and references therein. The theory is, however, plagued by many difficulties. 
One of them is that generical non-linear completions of the Fierz-Pauli linear theory propagate an 
additional scalar degree of freedom with various pathologies, the so called Boulware-Deser ghost [2]. 
Perturb at ively, the ghost manifests itself as a scalar field with a wrong sign kinetic term. Despite this, 
the ghost was argued [3] to play a role in the Vainshtein mechanism [3j of alleviating the so-called van 
Dam-Veltman-Zakharov discontinuity [5j. A related difficulty with massive gravity is that the theory 
becomes strongly coupled at a very low energy scales [6] . 

Recently de Rham, Gabadadze and Tolley [9] suggested that a certain non-linear completion 
of the Fierz-Pauli theory is free from the ghost degree of freedom. Moreover, it was argued that the 
Vainshtein mechanism is still at play, and that the cut-off is raised to a much larger energy scale, so that 
the ghost free massive gravity can in particular be trusted in the Vainshtein mechanism regional A non- 
perturbative formulation of the action was given in |11 4 I12 [ IT5 ] I14]. including its natural description as a 
sector of a bi-metric theory of gravity, with the action being the sum of two Einstein-Hilbert actions for 
the two metrics, g^ v and g~ v , plus a local interaction V(g +Pp g~ u ) between them. Linear combinations 
of the metric perturbations describe a massless and a massive spin-2 particle. Consequently, a number 
of papers have studied the proposed interaction [Ml HH [13 [IH HH [2QI EH [22l EH [24]. One of the 
key goals of these works was to establish the ghost-freeness non-perturbatively, by means of canonical 
analysis. However, the latter turns out to be rather cumbersome, both in the single massive gravity 
and in the bi- gravity frameworks, owing to the complicated expression of the potentials V argued to 
give a ghost-free theory. In particular, the potentials include square roots of metrics and did not seem 
to be natural from any geometric point of view. 

A new insight came from the work [25] (see also [26] ) where it was pointed out that these potentials 
take a much simpler form in the tetrad formalism. The tetrad is a set of four one-forms e 1 = e^dx M , 
related to the metric by 

= e i e tviJ- (!) 

The bi-gravity ghost-free interactions, rewritten in terms of the tetrads e+ and e{_ for the two metrics, 
reduce to the following 4-forms 

eijKLe T + A e J + A ef A e L _ , euKL^i A e J + A e5 A e L , eijKL^ Ae J _Ae*Ae L . (2) 

This is an important advance, which in particular shows that the ghost-free theory is not an artificial 
monstrosity, but a rather natural theory, once the correct set of variables is identified. The actual 
equivalence between this formulation and the original one in terms of metric variables is not exact, 
but depends on the vanishing of a certain antisymmetric contraction of the tetrads. This was shown in 
[25] to happen always at the linearized level, where the condition amounts to a gauge fixing. However, 
it is not true in general. Conditions for its validity were recently spelled out in [27] . 

Modulo the above discrepancy, the tetrad (or more generally vielbein) formulation also makes 
it completely transparent how the primary constraint responsible for removing the ghost degree of 
freedom arises. Indeed, as pointed out in [25], the interaction terms ([2]) are all linear in the lapse and 
shift functions of both metrics (due to the wedge product used in constructing them). Varying with 
respect to these lapse and shift functions, one gets the set of primary constraints. A part of this set is 
expected to be first class, and realize the algebra of diffeomorphisms of the bi-metric theory. Another 
part is expected to be second class. In particular, the constraints that are supposedly removing the 

1 See, however, a recent paper [10] for a reevaluation of the strong coupling issue. 
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ghost degree of freedom are a second class pair, with one of these constraints arising as primary and 
one as a conjugate secondary constraint. Unfortunately, the analysis quickly becomes rather messy, 
even in the simpler vielbein formulation of [25] • In the metric formalism, the (difficult) computation 
of the secondary constraint was carried out only recently in [T3] and in the Stiickelberg formulation 
in [28]. The first of these papers also gives an argument as to why the stabilization procedure of 
requiring the preservation of the secondary constraint by the time evolution just fixes one of the 
Lagrange multipliers and does not produce any new constraints. 

However, to the best of our knowledge, the complete Hamiltonian analysis has never been carried 
out in the first order formalism. Given that this is the formalism where the primary second-class 
constraints are immediate to see, one could possibly expect also some simplifications in the structure 
of the elusive secondary constraint. The aims of the present paper are to carry out such first order 
formulation analysis, for a model closely related to the one in [25], as well as to put the ghost-free 
massive gravity in the framework of the so-called Plebanski chiral formalism. 

Thus, in this paper we study a version of the ghost-free massive gravity using what can be called 
a "chiral" description, based on the chiral splitting of the local Lorentz gauge group. Such a chiral 
description is well-known in the case of pure gravity (with or without cosmological constant), and 
goes under the name of Plebanski formulation of general relativity [29, [30| 13 lj . It can be seen as a 
generalization of the tetrad formalism, where the fundamental fields are an SU (2) connection A and 
a two-form B valued in the su(2) chiral subalgebra of the Lorentz algebra. The SU(2) connection 
replaces the spin connection, and B replaces the tetrad one-form of the tetrad description. The metric 
is then constructed via the so-called Urbantke formula [32], 



where the epsilon symbol with an overtilde stands for the densitized anti-symmetric tensor, that does 
not need a metric for its definition. Much like the tetrad is a "square root" of the metric, the two-form 
here is its "cubic root". This apparent complication brings in several advantages. First, the phase 
space of general relativity can now be described in the same way as an SU{2) gauge theory, with 
the "triad" field for the spatial metric being canonically conjugate to the spatial part of the SU{2) 
gauge field. This simplifies the canonical analysis tremendously. Calculations can be performed very 
efficiently, and one is able to work out the complete constraint algebra with much less effort than 
in the metric (or vielbein) formulation. Secondly, it allows to automatically take into account the 
extension of the theory to a first order framework, where the connection is an independent field, like 
in the Palatini formulation of general relativity. 

The theory we consider in this paper is given by two copies of chiral Plebanski actions, coupled 
together in essentially the same way as was advocated in [25] in the context of the vielbein formulation. 
More precisely, the interaction we consider is induced by the term 5ijB +l A B~ 3 , where B + and B~ 
are the two-form fields for the two copies of Plebanski gravity. We study the canonical structure of the 
theory, and in particular, we explicitly compute the algebra of primary constraints and the resulting 
secondary constraint, and show that the system has seven degrees of freedom so that it is free of the 
Boulware-Deser ghost. We also prove that the stabilization procedure closes (under some assumptions 
to be clarified below). Thus, the results of this paper provide an independent proof of the absence 
of the ghost and moreover generalize it to the first order formalism. As an illustration of the power 
of our formalism, we use it to compute the Poisson brackets of two (complete, with all second-class 
constraints taken into account) smeared Hamiltonians. This is a useful exercise, as it computes the 
quantity that is known to have the interpretation of the physical spatial metric. As far as we are 
aware, such a computation has never been performed in any other formalism. 

The price to pay for the advantages of the chiral formulation is that, in the Lorentzian signature, one 
has to work with complex valued two-forms, and at some point supplement the theory with appropriate 
reality conditions. For the chiral description of pure general relativity these reality conditions are well- 
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understood. However, new issues arise after two copies of general relativity are coupled together, and 
for the interacting theory the problem of reality conditions has to be solved essentially from scratch. 
This is not attempted in the present work, see, however, some comments on this issue in the main 
text. In particular, an interesting observation is that the problem of finding suitable reality conditions 
seems to be related to the discrepancy between tetrad and metric formulations mentioned above. With 
this question lying open, our results and proofs apply to Euclidean signature bi-metric gravity. 

As will be shown below, the interaction term we consider corresponds to the second term in ([2]), 
whereas the first and third ones cannot be obtained in our formalism. On the other hand, the use 
of the 2-form chiral formalism opens a door to an infinite class of interactions leading to ghost-free 
bi-gravity theories, as opposed to the only three terms of ([2]). This comes from the possibility to 
modify the dynamics described by the Plebanski action without introducing extra degrees of freedom, 
see e.g. [33] for a reference. As we will argue at the end of the paper, such modifications can be 
incorporated in the coupled system and are also free from the Boulware-Deser ghost. While we will 
not study these theories here in any details, it is interesting to note their existence. 

The organization of this paper is as follows. In the next section we introduce the chiral Plebanski 
formulation of general relativity and demonstrate how simple is its canonical analysis. In section [3j 
we define the theory to be analyzed and discuss its relation to the metric formalism. Section [J] is 
the main part of the paper where we present the canonical analysis of the interacting system. The 
resulting canonical structure is summarized in subsection 14.51 whereas in subsection 14.61 we discuss 
an infinite parameter generalization which appears to preserve this structure. The final section [5] is 
devoted to conclusions. A number of appendices present some more technical results, and in particular 
the computation of the secondary constraint and the Poisson bracket of two smeared Hamiltonians. 



2 Summary of the chiral Plebanski formulation 

2.1 Plebanski formulation of general relativity 

The chiral Plebanski formulation of general relativity is described by the following action [29] 

1 



S[B,A,9] 



B A F 



W3 - \5V)B l AB : > 



(4) 



The indicedH i,j = 1,2,3 are suc(2) ~ so(3,C) ~ sl(2,C) Lie algebra ones, with the Lie algebra 
viewed as a vector space of complex dimension 3. The field A 1 is the connection one form, F % = 
&A l + £ l i k Ai A A k is its curvature, and the field B l is a Lie algebra-valued two-form. A is a multiple 
of the cosmological constant, and the Lagrange multiplier field \E'* J is required to be traceless. Its 
variation enforces the so-called metricity (or simplicity) constraints 



B i AB j = - 5 ij B k A B k . 
3 

A solution to these constraints can be conveniently written using tetrad one-forms valued in so (4, 
as 

3±0i (J A J 



(5) 



B l = P, 



(e 1 A e J ) , 



(6) 



where I = 0, i and Pt z j KL are two chiral projectors defined in terms of the internal Hodge star * as 



I J : =o( 1±CT *)/J = 7 \°I 6 J ~ 6 I 6 J ± <rej j ) 



(7) 



2 Since the su(2) indices are raised and lowered with the unit metric, we do not strictly follow the rule that the indices 
one sums over should be in opposite positions. All repeated indices are assumed to be summed. 
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Here a = i for the Lorentzian signature, and a = 1 for Euclidean. Plugging the solution ([6]) back into 
the action Q, one obtains the Einstein-Cartan action plus a boundary term. The equivalence with 
the Einstein-Hilbert action then follows as in the Palatini formalism: one varies with respect to the 
independent connection, obtaining an equation that identifies it in terms of the unique torsion-free 
Levi-Civita connection. 

The action given above uses complex fields and therefore describes complexified general relativity. 
To get a real theory, one has to impose appropriate reality conditions. These are easiest to state by 
saying that the metric ([3]) should be a real metric of the required signature. Such conditions can be 
given directly in terms of the fields that appear in the Lagrangian. For metrics of Euclidean signature 
the reality condition is simply that all fields are real, and thus the connection and the B field are 
so (3, R)-valued. For Lorentzian signature the conditions are quadratic in the fields, and can be stated 
as 

B i A {B j f = 0, Re(B i A B' 1 ) = 0. (8) 

First, one requires that the complex conjugate of any of the 3 B J, s is wedge-orthogonal to any other 
B l (or to itself) which gives nine real equations. In addition, one demands the 4-form obtained as the 
trace of the wedge product of the -B's be purely imaginary. The reason for these conditions is clear 
from ([6]): when the B i are constructed as the self-dual part of the two- form e 1 A e J , with e 1 a real 
Lorentzian signature tetrad, its complex conjugate is the anti-self-dual part, and both conditions in 
© holdl 

The reality conditions for B l also induce reality conditions on the connection A 1 . These are 
obtained from the field equation for the connection that follows from (jl]). The field equation can be 
solved for A 1 in terms of B l , and then the reality condition (|8j) implies a reality condition on the 
connection. Alternatively, one can state the reality condition on A 1 by passing to the Hamiltonian 
formulation, as shown below. 



2.2 Canonical analysis 

As a preparation for a more complicated canonical analysis of a bi-gravity model, we recall here the 
canonical analysis of the chiral Plebanski formulation @. It starts with the usual 3 + 1 decomposition 
of the action which leads to the following Lagrangian 

C = \ ~e abc (Bi a Fl + B\ c Fl a - (¥* - A^')«) . (9) 

It contains only one term with time derivatives, thus the phase space of the theory is parametrized by 
the spatial components of the connection A l a , with conjugate momentuir@ 

~ f)C 1 

The remaining fields Aq, B^,^^ are non-dynamical. Since they enter the Lagrangian ([9]) linearly, 
they are Lagrange multipliers generating constraints. To simplify the analysis, it is convenient to 
decompose B l 0a as 

Bl a ={N5^ + e^ + lp)El (11) 

where N,U,U 3 describe its trace, antisymmetric, and tracefree (b n = 0) parts, respectively, and we 
assumed the invertibility of E ai denoting by 



K = \e^ k E)Ei (12) 



3 When describing the Lorentzian theory, it is also convenient to put an imaginary unit factor i in front of the action, 
as to get precisely the Einstein-Cartan action in the end. 

4 The tilde above the symbol keeps track of the density weight of this pseudotensor. We will also use a tilde under 
the symbol to characterized negative weight pseudotensors. 
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its (densitized) inverse. Invertibility of the triad is a necessary condition to reproduce general relativity, 
perfectly analogous to the invertibility of the metric. After the decomposition (|lip . the constraints 
obtained by varying the action with respect to $? tJ simply imply 

Ifi = 0. (13) 

Thus, these components of the Bq g field can be omitted from now on. In terms of the remaining fields 
and using the notation 

N a := VE?, (14) 

the above Lagrangian can be rewritten as 

C = Efd A\ + J%D a E* + N a E b F* ab + - e» k [E^F^ + Xe&EftfEty , (15) 

where D a is the covariant derivative defined by the connection as D a X l = d a X l -\-e^ k A J a X k , and we 
omitted a total derivative term. This representation identifies Ef as the densitized triad determining 
the spatial metric, with the last term in the Hamiltonian being the cosmological term, cubic in the 
triad and thus giving the volume form of the metric. Notice also that the covariant expression (|6|) is 
recovered from ()10|) and (|lip upon using the constraint fjl3[> . definition (|14H and the standard ADM 
decomposition of tetrad. In particular, N and iV a are identified as the usual lapse and shift functions, 
respectively. 

The resulting Lagrangian identifies the symplectic structure to be given by the Poisson bracket 

{A\{x),E b {y)} = 5 b a 5)~5 3 {x-y). (16) 

This is the same phase space of an SU(2) gauge theory, which is the main advantage of this formulation, 
and also the basis for Ashtekar's approach and loop quantum gravity [341 ES]- The constraints arising 
by varying with respect to Aq,N, N a are the Gauss, Hamiltonian and diffeomorphism constraints, 
correspondingly. They are first class since they form a closed algebra under (I16p . To write the algebra 
explicitly, we define the smeared constraints 

Q<t> ■= J d 3 x0 l D a E°, 

V n := J d 3 x N a E u F l ab - G N a Aa = J d 3 x N a (E bi d a Al - d b (E bi A^ , (17) 

U N := \j d 3 xNe^ k (EtE b F k b + Xe^EfE^) . 

Notice that we have shifted the initial constraint coming from N a by a term proportional to the Gauss 
constraint. This is convenient because the resulting constraint precisely coincides with the generator 
of spatial diffeomorphisms. Then it is straightforward to verify the following relations 

{GfajGfa} = Gfaxfa, G<(,} = GN a d a fa {G^Un} = 0, 



where we used notations 



(a x by :=e ijk aib, 
L a (N, M) 



L(N,N) 
N a (N, M) 



•j °fc i 

N b d b M a - M b d b N a 



(19) 

N a d a N - Nd a N a , 
EfE b (Nd b M - Md b N) . 
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The above algebrsd ensures that all constraints are preserved under the evolution generated by the 
Hamiltonian 

Htot = —GA +N a Aa — ^a? ~~ H- N > (^) 
which is given by a linear combination of first class constraints, as it should be in any diffeomorphism 
invariant theory. 

As a result, the constraints generate 3 + 3 + 1 = 7 gauge symmetries which reduce 9 + 9 = 18 
dimensional phase space to a 4-dimensional one which describes two degrees of freedom of a massless 
graviton. It is worth stressing how much simpler is the above Hamiltonian analysis as compared to 
the usual metric ADM analysis, or to an analogous analysis using tetrads. It is this simplicity that we 
would like to take advantage of, and write a similar "crural" model of the ghost-free bi-metric gravity. 

Finally, we briefly discuss the reality conditions in the canonical picture. In the case of the 
Euclidean signature, everything is real. For the Lorentzian case, the spatial triad Ef is taken to 
be real, which ensures the reality of the spatial metric and can be seen to be a part of the reality 
condition ([8]) discussed above. The remainder of the condition ([8]) becomes the requirement that the 
lapse N is purely imaginary and the shift N a is real. The reality condition for the connection can 
then be obtained by requiring that the reality of Ef is preserved by the time evolution (generated 
by the constraints). This fixes A l a to be of the schematic form A l a = T l a -\- \K l a , where r„ is the 
connection compatible with the triad, and K\ is the extrinsic curvature of the spatial slice embedded 
in 4-geometry. For more details on this, the reader can consult |34j . 

3 A chiral bi-gravity model 

A bi-gravity theory is a model represented by two copies of general relativity coupled by a certain 
interaction term. Our idea is to use the chiral description of general relativity given by the Plebanski 
action (j4j). The simplest interaction then leads to the following action, 



S[B ± ,A ± ,^ ± ] = J 



B +i A p+i + B -i A p-A 

(21) 

~ \ + S ij )B +i A B +j - - \-5 ij )B- 1 A B- j + 2aB +i A B~ l , 

where the indices ± distinguish the fields from the two sectors. The coupling constant a will later on 
get related to the mass of the second graviton. 

In absence of the interaction term, i.e. for a = 0, the gauge group of the theory is given by the direct 
product of the symmetries of each of the copies of Plebanski theory, which include diffeomorphisms 
and local gauge rotations. The interaction breaks the total group Diff + x SO + (3) x Diff _ x SO_(3) to 
its diagonal subgroup, and it is this fact that is responsible for a larger number of propagating modes 
in (|2ip as compared to the non-interacting theory. 

Our model can easily be reformulated in terms of the tetrads instead of the two- forms B^. Indeed, 
in each ±-sector the variation with respect to the Lagrange multipliers <I' ±l - ? imposes the simplicity 
constraints ©, which imply that B^ are the self-dual projections^ of the two-forms e± A e± as in 
([6]). Thus, each Plebanski term reduces to the Einstein-Cartan action for the tetrad with its own 
cosmological term. At the same time, the interaction term takes the following form 

77 euKLe{ A e J + A e K A e L + e+ A e{ A e_j A e_j. (22) 



5 The constraint algebra offers a direct way of confirming the triad interpretation of Ef. Indeed, as nicely explained in 
5], the spatial metric can be read off the right-hand-side of the Poisson bracket of two Hamiltonian constraints. Namely, 



it should appear in the function multiplying the diffeomorphism constraint. From the last line of (fj9j), we observe that 
this quantity is given by EfE", which makes clear that the spatial metric is constructed from Ef as a triad. 

6 In © the indices ± denote the chirality and should not be confused with the indices ± of this section which 
distinguish the two Plebanski sectors. 
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The first contribution is exactly the symmetric interaction term in ([2]) which has been argued to 
generate a ghost-free massive gravity [25]. This is why we take (|2ip as our starting point. Notice also 
the presence of the second contribution, not of the type ([2]) considered in [25J. It is consistent with all 
symmetries of the theory, except it is parity-odd, unlike the rest of the Lagrangian. 

On the other hand, the relation of our model to the metric formulation is more subtle. The reason 
for this is that the equivalence between the polynomials (|2|), and the ghost-free interactions between 
two metrics originally proposed in [IT] and related references, holds only if the tetrads satisfy the 
symmetry property 

J+^-vM = 0. (23) 

In [25] it was shown that this condition can always be realized in perturbation theory, where it basically 
amounts to a choice of gauge. However, as it was recently pointed out in [27], it does not hold true 
in general. Hence, already the tetrad formulation turns out to be slightly different from the metric 
formulation. Our model introduces a further difference due to the presence of the second contribution 
in (I22p . In our formulation, this term simply can not be discarded, as the whole structure (I22p 
comes at once. It is interesting to note that this contribution vanishes provided the condition (|23p 
holds. This fact puts our formulation of the ghost-free bi-gravity theory on the same footing as the 
tetrad formulation as both of them coincide with the one originally proposed in [71 [11] only for the 
configurations satisfying (I23p . 

Let us also make a few comments about reality conditions. In the case of Euclidean signature, as 
in the usual Plebanski theory, all fields can be taken to be real and one finds real Euclidean bi-metric 
gravity. In contrast, in Lorentzian signature, the situation with the reality conditions in the presence 
of the interaction term becomes more intricate than in the single gravity case. The simplest idea is to 
impose the same reality conditions which have been discussed in section I2TT1 in each of the two sectors. 
This renders each metric real and ensures also that each of the Plebanski actions is real modulo a 
surface term |37] provided one multiplies the whole action by a factor of i. As for the interaction 
term (|22p . the extra i cancels the factor of a = i in the first term in (|22p and makes it real as well. 
However, the second term in (]22p then becomes purely imaginary, thus the total action is not real. 
This implies that for Lorentzian signature the standard reality conditions of Plebanski theory are not 
appropriate and should be modified to deal with the bi-gravity case. While it is possible that there 
is some more sophisticated choice, we do not consider this issue any further in the present paper, 
and content ourselves to perform the analysis for the Euclidean signature when all the fields are real 
(or, equivalently, for the complexified theory, when all fields are complex). Nevertheless, let us note 
that, as pointed out in the previous paragraph, the term spoiling the reality of the action vanishes on 
configurations satisfying (|23f) . As a result, we come to an interesting observation that the condition of 
the reality of the action (for the standard reality conditions on the fields) precisely coincides with the 
condition for the equivalence with the metric formulation. However, just like it is not obvious how to 
include that restriction from the start in the tetrad formalism, it is also not obvious how to achieve 
the reality of pip in a simple geometric way. 

4 Canonical analysis and ghost-freeness 

In this section we perform a careful canonical analysis of the model (I2ip . obtain explicit expressions for 
all constraints, including the secondary constraint responsible for the absence of the Boulware-Deser 
ghost, and compute their algebra. The results of this analysis are summarized below in subsection 

EE 
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4.1 3 + 1 decomposition and field redefinition 

The Hamiltonian analysis of the action (|21|) can be done in the same way as in section 12. 21 The first 
step is to perform the 3 + 1 decomposition which gives rise to two actions of the form for the 
±-fields, plus the term coupling the two sectors 

Ant = a (BgB£ + BgBg) ~e abc . (24) 

In both sectors we can introduce the fields E, N and N a as before in (IIOD and (llip . The two traceless 
components are again constrained to vanish by the simplicity constraints generated by ^ ,±JJ . In 
terms of the remaining components, the interaction term (I24p takes the following form 



C int = a ( N+ea bc ^ k E+ b E+ c + 2ea bc N +b E +ci ) E~ a ^ 



3 



+ a ( N-e^ E~ b E^ + 2^"^ ) 



The canonical form of the action is then the sum of two copies of (|15j) plus the above interaction. 
Being linear in lapse and shift, the interaction term directly contributes to the constraints already 
present in the ±-sectors. As will be shown below, these contributions break the "off-diagonal" part of 
the gauge symmetry of the non-interacting theory, leaving only its "diagonal" part as the symmetry 
of the full theory. For this reason it is desirable to introduce variables that will be adapted to this 
pattern. Motivated by this, we do the following field redefinition: 

pa J?+ a _L_ /? — a F+ a F — a 

h i h i + h i > H i &i ~ &i > (r,(! a \ 

Ai = {A? + A?)/2, vt ■= (Af-A~>)/2. ^ } 

Of course, only the variable A l a remains a connection, with the field rf a transforming as a matter field 
under the diagonal SO(3). The curvatures F ± decompose as 

where D a is the covariant derivative with respect to the connection A a . Performing a similar redefi- 
nition for the Lagrange multipliers 



N a 
N 



(A+* + v)A r 

(N +a + N- a )/2, G a 
{N + + N-)/4, G 



(Af - V) /2, 

(N +a - N- a ) /2, (26b) 
(N+-N-)/A, 



allows to disentangle the diagonal and off-diagonal sectors, thereby making the diagonal symmetries 
manifest. Distinguishing the diagonal and off-diagonal constraints by putting a hat on the latter, the 
total action after the change of variables takes the following canonical form 

S = Jdt d 3 x {e^OqA^ + H ai d Q rf a + 4> l Qi + ^Qi + N a C a + G a C a + NU + G%) . (27) 

Before we give the explicit form of the constraints, it is useful to contrast the above change of 
variables with what is usually done in the metric formalism. There one keeps the original metrics as 
fundamental variables, and mixes only the Lagrange multipliers (using "geometric averages" rather 
than algebraic ones as above), see e.g. [55] . However, upon linearization the mass eigenstates turn out 
to be the sums and differences of the (perturbations of the) two initial metrics, h£ v ± h~ v . Extending 
this redefinition to the full theory, one may work with 

9/j.l* = g^,v + 9fiui ifj,v = g^u - g^v ( 28 ) 
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This is morally what we are doing here. In particular, we expect that perturbations of Ef describe the 
degrees of freedom of a massless graviton, and those of Hf a massive spin 2 field. However, there is 
no reason to expect that the linear transformation (|26p or (|28|) decouples the physical metric and the 
massive field beyond the linearized theory. In fact, even without discussing the coupling with matter, 
the physical metric can be identified by computing the Poisson bracket of two smeared Hamiltonian 
constraints, as explained in footnote [5j For our model this calculation is performed in appendix ICl and 
shows that the physical spatial metric is a non-trivial function in phase space, and does not coincide 
neither with g ab defined by Ef, nor with the ones defined in (|28p . The same situation arises in the 
metric formalism, see e.g. |39j . Thus, disentangling the relation between the physical spatial metric 
(as appears in the algebra of diffeomorphisms) and the phase space variables is a complicated task, 
whatever formulation is used. Despite the redefinition (|26p does not fulfil this task, it does disentangle 
the constraint algebra and therefore is very convenient. 

Let us now discuss the two Hamiltonian constraints which read as 



U = l e ijk (E ai E bj + H ai H bj )(F k b + e h 



Mm rf a vT) + e ljk E ai H bj (D a rj k 



DbV k a ) 



a + ^j EfE b E% + 3/3E?EjH%. + 



a 



fpa rrb ttc i ft rja rrb ttc 
&i H j H k + 2 H i H j H k 



(29) 



and 



H = e ijk E ai H bj (F k b + e 



klml „m\ 



£1 E-E)Et + la + - EtE)H% + -f E?H*H c k 



- - a I HfH b H c k 



(30) 



where we denoted 



A 4(A+ + A-), /3 = I(A+-A-). 



(31) 



2 v ' 2 

These expressions, and the subsequent analysis, can be significantly simplified by restricting the pa- 
rameters as 

A = -2a, = 0. (32) 

These restrictions can be justified as follows. First, the combination a + A/2 can be identified with 
the effective cosmological constant of the combined system. Since the presence of the cosmological 
constant has a very little effect on the constraint algebra, we can safely set it to zero. The second 
restriction arises from the observation that it removes from the Hamiltonian constraint % the term 
linear in H whose presence would imply that the Lagrangian for the massive field contains a tadpole 
and is not in its canonical form. 

After we have made the choice (132p . one remains with only one free parameter a, which can be 
identified with the mass of the H field. The full set of simplified constraints resulting from this choice 
can be found in (l34l) and will be discussed in the next subsection. 



4.2 Primary constraints 

The canonical form of the action (|27p shows that the phase space of the theory is spanned by the 
fields A, E, ij, H and carries the (pre-)symplectic structure encoded in the following canonical Poisson 
brackets 

{Ai(x),E b (y)} =S b 6yS 3 (x-y), {rf a (x\H b (y)} = ^5 3 (x - y). (33) 
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The other variables introduced in (|26b|) are the Lagrange multipliers for the primary constraints which 
have the following (smeared) expressions 



d 3 x^ 



d 3 x^ 



d 3 xiV a 
d 3 xG a 
d 3 xN 



D a E? + e ljk rf a H ak 



D a H? + e ijk if a E ak 

El(Fi b + e^ k rf aV k b ) + fl*(D a i£ - D b rf a 
Et{D a r,l - D b ri) + + ^Va^) - 2ae abc £f #f 

\ e^(EtE) + i^)(F a fc b + e Wm f/ a C) + e^EfH^D^ - D b r, k a ) 



(34) 



d 3 xG 



In complete analogy with the single copy of Plebanski theory, it is convenient to shift the constraints 
Cff and by a linear combination of the two Gauss constraints. Namely, we define 



e ijkE t H b AF k b + e klm rj l aV r) + e^{EfE] + H?ffi)D a rfi 



V 



Cfi — GN a A a — Gn 



d s xN a 



E\d a A\ - d b ( E\A\ ) + H?d a 4 - d b ( H°nl 



N 

:= - g N a Va - N a Aa = I d 3 x N a E\d a rf b - d b ( E\rf a ) + H\d a A\ - d b ( H\A\ ) (35) 



2a / ci 3 xiV a 



%s>i>cE h iHf 



It is now a straightforward although tedious exercise to compute the algebra of the primary con- 
straints. In contrast to the case of pure gravity, it is not anymore first class: some of the commutators 
acquire contributions which are non-vanishing on the constraint surface. The first class part of the 
algebra is as follows, 



{G<f>,T-tN} 



^ix02i {G,p,G^} = G<f>xip, 

GN°*da<t>-> {DffiG^} = GN a d a tl>-> 

0, {G^Hg} = 0, 



n 



L(N,N)> 



{d^Ug} = n 



L(N,G)> 



{GipnGfa} — Gfaxfai 
\Pg-,G<]?$ = Gc a d a <i>, 



L(N ,M) ' 



(36a) 



L(N,G) ' 



and 



{Hn,7-Lm} - ^v(N,M) + ^U(N,M) + Gv a (N,M)A a +U a (N,M)Va + Gv°>(N,M)Va+U«(N,M)A a > 
{H-GiUf} = ^viCF) + T> U(G,F) + < 3v a (G,F)A a +Ui(G,F)r) a + Gv a (G,F) Va +U a (G,F)A a , 

where we used the notations introduced in (|19p as well as 

V a (N, M) := (e?E$ + HfH^j {Nd b M - Md b N) , 
U a (N, M) := (E?Hl + HfEf) (Nd b M - Md b N) . 



(36b) 
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The constraints Gif> and T>^ weakly commute with all other constraints and form exactly the same 
subalgebra as in (fTSj) . It is thus clear that they represent the generators of the usual gauge and 
diffeomorphism transformations. 

The remaining Poisson brackets are given by 



{v d ,n N } 
{n N ,u G } 

where 



G G °8 a i> + 2a J (Px^G^e^ (E)El - HjHty , 



-4a / iPxNiFyJ^E^H?, 
-4a / SxG^eabcHfE^H^ 



V l { g,f) + 4a / d 3 xe flbc G a F b (Efd d Et - Hfd d H?) , 



^L(G,N) + a ^i(Af,G) a ^ch(N,G) 



(36c) 



+ 4a / d 6 xNe, 



G d A a EfE b j H] + G a Vd [E% [E]E<f - HfHf) + HfEjHj 
aQ 



?C T?d 



TC Tjd 



jd i?b TTC 



^L(G,G) + a G<f> 2 (G,G) 



+ 4a / (Pi&pfc 



Vi(G.G) 
ra fib Tie 



G a A l d H a i E)B c j + G a rf d [ H? lEfE* - H°H?) + EfEjH 



?d fpb : 



~ V U{N,G) ~ V V(N,G) ~ Gu"(N,G)A a +V"(N,G) Va 



Qua 



U a (N,G) Va +V a (N,G)A a 



2a / d*xNGe abc eV k (EfE^E^Ef + H?H°H c k Ht - EfH^Ef - HfE°,E^Hf ) rf d 



b 7jc ; 



pa zjb tjc r?d 



ra r?b i 



b\(N,G) :=2NG a e i > k e slbc E)Hl 

%(N, G) := NG a e^ c + H)H% 



(37) 



As is expected, all contributions non-vanishing at the constraint surface are proportional to the mass 
parameter a. 



4.3 Secondary constraint 

Since the primary constraints do not form a closed algebra, the Dirac's stabilization procedure does not 
stop at the first step and we have to study the conditions ensuring that the time evolution preserves 
the constraints. The evolution is generated by the total Hamiltonian given by a linear combination of 
the primary constraints 

fltot = -Q$ -G^-V^-V 3 -U N - H G , (38) 

where we denoted = + N a A a + G a rj a and ip = ip + N°"q a + G a A a . Since G% and T> a weakly commute 
with all primary constraints, the requirement of their stability under evolution does not generate any 
conditions. For other constraints one finds 

ti ~ {f> 3 , d} + {H N , Gi} + {n G , Gi] « o, 
v a ~{G^A} + {f>G,Va} + {n N ,A} + {n G ,v a }~o, 

(39) 

H*{gj,H} + {V d , H} + {H G , H} « 0, 
ii ra r, H] + {V d , H] + {U N , H} ra 0. 
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The first two equations can be used to fix the Lagrange multipliers for the "off-diagonal" Gauss and 
diffeomorphism constraints^ 



G a 

ft 



N{n,g l } + q{n,g t }) 

n ^{n^jUGjM^iAM + {n,v a }) {A*, ft}" 1 
+G({n,g j }{g J ,v b Y 1 {v b ,v a } + {n,v a }) 



(40) 



Plugging these expressions into the third equation, one observes that the terms proportional to N 
cancel and the whole expression is proportional to G. As a result, it leads to the following secondary 
constraint 



V = {H, H} + {H, Gi}{gi,A}~ l {A, H} + {% Va}{V a ,Qi}~ l {Qi,n} 



(41) 



On the surface of this constraint, the last equation in (|39[) is satisfied automatically. Thus, is the 
only secondary constraint arising at this stage of the Dirac's procedure. 

Having calculated in ()36c[) all non- vanishing Poisson brackets entering the expression for ^f, one 
can proceed evaluating this constraint explicitly. This is done in appendix[A]where the following result 
is obtained: 



-4a 



e[Ef- HjEfoHi) +h(Hf- E*H{E b ) rf a + A^T 



where 



det-E?, 



A 

/; = 



= (e - hT)e -(h- eT)h, 
detm, T 



Tji jpa 



rpi zra 
E a H i > 



(42) 



(43) 



and 
T 



e I E%E b k + H a k H b ) +(h- eT)H a k E b e&EfHl&grt [Hf (4 E] - Hj H r A + E\E]E\ 



h(E a k E b + H a k H b 



+ (e - hT)H a k E b e^EfHfe^fvt. 



Efm^ 



+Cabre 13 Vl ( E i H k ~ H i E k) tgcd { E k' E k' + H k' H k>) E l H l &pq H m E L E n H n- 



(44) 



Although the result looks very complicated, the constraint features the simple properties of being linear 
in the field rj, and not containing any derivatives nor dependence on the connection A. The expression 
(|44p can be further manipulated using various identities and eliminating the epsilon tensors, but the 
result contains a number of monomials and is not particularly enlightening. Finally, it is interesting to 
note that T is proportional to the combination e^^E^Hf, analogous to the antisymmetric combination 
of the tetrads (j23|) which must vanish for having agreement with the metric formulation. Therefore, 
on configurations satisfying this condition, the secondary constraint crucially simplifies. 



4.4 Stability condition for the secondary constraint 

The next step is to study the stability of the secondary constraint ^> . However, due to its complicated 
expression, it is difficult to do this by a direct computation. Nevertheless, some important conclusions 
can be made if one realizes that \P is given by a "partial Dirac bracket" of the two Hamiltonian 
constraints. By partial Dirac bracket we mean here the Dirac bracket constructed using the minimal 

7 In the following equation, the Poisson brackets with a slight abuse of notation denote the expressions appearing 
under the integral in the part of the commutator non-vanishing on the constraint surface. In other words, we consider 
the Poisson brackets of the non-smeared constraints dropping the distributional 5-factor. 
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set of non-commuting primary constraints. Since generically this is just a subset of the full set of 
second class constraints of the theory, it is not yet the final Dirac bracket. Nevertheless, it is a useful 
object since the final Dirac bracket is obtained from it using the missing non-commuting constraints 
in the same way as it is constructed from Poisson bracket using all second class constraints. 

In our case, the minimal set of non-commuting primary constraints consists of Q and T>. Their 
Dirac matrix and its inverse are given by 

{A,Gj} {A, A})' V -{AiGjY 1 o )■ 

The partial Dirac bracket is the standard Dirac's formula 

{A, B}' D = {A, B} - {A, £ S }W S '{<£ S , , B}, (45) 

where <t s = (Gi,T> a ). It is immediate to see that the expression (|41[) for the secondary constraint is 
equivalent to 

VP = {%x,%Y D) (46) 

where we smeared the Hamiltonian constraint T~L with the trivial smearing function 1 just to remove 
the (5-function factor from the right-hand-side. Besides, it is useful to note that the total Hamiltonian, 
after plugging in the expressions (|40p for the Lagrange multipliers conjugate to the constraints £ s , can 
be written as 

H t ot = -g^-T> g -n' N -n' G , (47) 

where W, n are the two Hamiltonian constraints corrected by £ s such that they weakly Poisson 
commute with all other primary constraints. Due to this, the stability condition for \1/ takes the 
following form 

y = {H N ,yy D + {HG,y}'D~0, (48) 

where we used the fact that ^ weakly commutes with Qi and T> a as well as that under the partial 
Dirac bracket Ti' and n can be replaced by the original constraints. 

The stability procedure ends if (|48p can be interpreted as an equation fixing one of the Lagrange 
multipliers. A necessary condition for this is that it should be algebraic in N and G, i.e. it should 
not contain their spatial derivatives. It turns out that this property follows just from the fact that ^ 
is well defined, i.e. that {Hn^gVd does not contain derivatives of the smearing functions. Indeed, 
let us write the most general form of the commutator 

{n N ,{ni,ny D y D = x a d a N + Ny, (49) 

where X a and y are functions on the phase space. Then, using Jacobi identity, one obtains 

{u N , *y D = {Hi, {M N ,ny D y D + {h, {H u n N y D y D « n {x a d a \ + y) = Ny, (50) 

where at the second step we used (|49p in the first term and the fact that any Dirac bracket with a second 
class constraint vanishes to remove the second contribution. A similar computation can be done for 
the second term in (f48l) if one notice that (l46j) can be equivalently rewritten as $ = {H,Hi}' D . Thus, 
we proved that the stability condition (|48p for the secondary constraint does not contain derivatives 
of N and G, and it is of the form 

Ny + qy = o. (51) 

For generic configurations in phase space, the functions y, y are non-vanishing, hence the stability 
condition (|5ip of the secondary constraint is a linear algebraic relation between the Lagrange multipli- 
ers iV and G. As a result, the Dirac's procedure stops at this point, there are no tertiary constraints, 
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and the secondary constraint \& is second class conjugate to one of the two Hamiltonian constraints. 
This provides the complete description of the canonical structure of our model in generic case. 

One may consider however the possibility of having configurations in phase space for which the 
two factors in (|51j) vanish simultaneously on the constraint surface, 

y = {^i,*}ij«o, y = {Hi,*}' D * o. (52) 

If this is the case, ^ becomes weakly commuting with other constraints and thus first class, further 
reducing the number of degrees of freedom. Although for a generic phase space configuration this does 
not happen, there still can exist some subsectors where ()52p holds. In fact, in appendix [B] we show 
that such a possibility is indeed realized and there are common solutions to the total set of constraints 
and the conditions (|52p . In particular, we study a perturbative expansion around H a = 0, evaluate 
explicitly the secondary constraint in the quadratic approximation, and the stability condition to the 
linear order in H a . A common solution is shown to exist in the sector with H a = and the connection 
A a restricted to have a constant curvature set by the mass parameter a. However, the subsectors thus 
found appear to be not stable under the evolution, which makes impossible to interpret the conditions 
(|52p as consistent additional constraints. It is likely that they do not have a physical significance and 
arise due to our incomplete understanding of the decoupling of the diagonal and off-diagonal degrees 
of freedom at the non-linear level. 

4.5 Summary 

Let us summarize what we have found so far. Generically the theory possesses 15 constraints. Among 
them there are 7 first class constraints: Qi, T> a and T~L' N + T-L' G , with the Lagrange multipliers N and G 
related by the condition (|5ip . These constraints generate the local symmetries of the chiral Plebanski 
formulation of general relativity: SU(2) gauge transformations, spatial and time diffeomorphisms. 
The remaining 8 constraints are of second class: Qi, T> a , and either % or % depending which of 
them does not commute with vl/Jf| Altogether these constraints fix 2 x 7 + 8 = 22 variables in the 
initial 4 x 9 = 36 dimensional phase space. As a result, the theory turns out to describe 7 degrees of 
freedom, and it is free from the Boulware-Deser ghost. 

Expanding around a bi-flat background, the degrees of freedom can be identified with a massless 
and a massive spin-2 fields. Specifically, the massless graviton is encoded in the "diagonal" (E, A) 
sector supplemented by 7 first class constraints, whereas the 5 degrees of freedom of the massive one 
are carried by the traceless symmetric part of the "off-diagonal" (H, rj) sector, with the antisymmetric 
and trace parts being fixed by the second class constraints (see appendix [B]). As anticipated, the 
identification of the physical fields at the non-linear level is more complicated. An explicit calculation 
in appendix C shows that the physical spatial metric identified by the Poisson bracket between two 
Hamiltonians is quite a non-trivial function lacking a simple expression in terms of our basic variables. 

It is interesting to compare the resulting canonical structure of this model with the one of bi-metric 
gravity in the metric formalism. In the latter case, one starts with two uncoupled copies of gravity, 
with the group of gauge symmetries being two copies of the group of diffeomorphisms Diff, and as a 
result 2+2 propagating degrees of freedom. After the two copies of gravity are coupled, the group of 
gauge symmetries is broken down to the diagonal group of diffeomoprhisms. This removes 4 of the 
gauge symmetries of the original uncoupled theory, and thus adds 4 propagating degrees of freedom. 
If, however, the interaction term is tuned to be of a special form such that it is linear in the two 
lapse functions, in addition one generates a couple of second class constraints (one primary and one 
secondary) removing one degree of freedom corresponding to the Boulware-Deser ghost and leaving 
us with 7 degrees of freedom. 

8 If they both do not commute with the secondary constraint, any of them can be chosen as second class. This is 
because they can be related to each other by adding a first class constraint, which represents a general ambiguity in the 
Dirac's approach. 
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A similar mechanism is at play in our model as well, but with some important differences. Like 
in the metric formulation, the group of gauge symmetries of the two uncoupled copies of Plebanski 
theory is also broken down to the diagonal subgroup by the interaction term. However, as (j'25p 
shows, the interaction term is linear in both, lapse and shift functions. Thus, while the off-diagonal 
group of gauge transformations is broken, there are still primary constraints corresponding to these 
transformations. The primary constraints for off-diagonal spatial diffeomorphisms together with the 
primary constraints for off-diagonal SO(3) rotations then turn out to form a second-class pair. Thus, 
even though by breaking the spatial diffeomorphisms and SO (3) rotations we have removed 3+3 gauge 
symmetries, we have only added 3 degrees of freedom. To say it differently, the interaction term just 
converted the corresponding constraints from first class to second class. As a result, they now remove 
a "half of the configurational degree of freedom each, as compared to removing one when they are 
first class. On the other hand, the removal of the Boulware-Deser ghost proceeds in the way analogous 
to the metric formulation and is based on the existence of the additional second class pair. While very 
similar picture was outlined already in [25], in this paper we have given its explicit realization for the 
first time. 

Finally, we have also observed that there are sectors in phase space, characterized by vanishing of 
the two (partial Dirac) brackets (|52p . where this classification of constraints may fail, in particular, 
leading to fewer degrees of freedom. In fact, the existence of special sectors in phase space with a 
drastically different canonical structure is not an unusual situation. It happens, for example, in the 
first order tetrad formulation of general relativity, and as well in the Plebanski formulation, where 
one can allow degenerate tetrads with the vanishing determinant of the metric. For instance, the 
degenerate sector of the non-chiral Plebanski formulation corresponds to a topological theory with no 
local degrees of freedom [40] . In principle, something similar might happen in our case as well| but 
the geometric interpretation of such subsectors is far from clear. In any case, they can lead only to 
a reduction of degrees of freedom, and therefore our analysis ensures that the scalar ghost has been 
removed throughout all the phase space. 

4.6 Modified Plebanski theory 

A striking property of the Plebanski formulation of general relativity is that it allows an extension 
to an infinite class of actions all propagating only 2 degrees of freedom. This should be compared 
with traditional modifications of general relativity like higher-order f{R) theories, which typically 
have additional degrees of freedom. Such special modification is achieved by turning the cosmological 
constant in (j4j) into an arbitrary functional A(\P) of the Lagrange multiplier The resulting action 
describes a modified theory of gravity propagating two degrees of freedom only [UJ . It can be argued 
to be related to an infinite sum of higher order curvature terms, organized as to not propagate extra 
degrees of freedom [12]. The specific form of the infinite summation depends on the specific choice of 
A(*). 

Relying on these results, one can suggest a similar generalization of our model (|2ip where the two 
cosmological constants are replaced by arbitrary functionals A =t ( l I /: ' = ). The argument implying absence 
of extra degrees of freedom goes through for both sectors of the coupled system, and the interaction 
does not appear to spoil it. Moreover, the interaction term is still linear in both lapses and shifts, 
which indicates the presence of the constraints responsible for the absence of the Boulware-Deser 
ghost. Therefore, even without performing a detailed canonical analysis, we argue that there is in fact 
an infinite class of ghost-free bi-metric gravities, characterized by the functionals A =t (^ /:l= ). This might 
be particularly useful concerning phenomenological applications of the theory. We leave the study of 
these issues for future research. 

9 Similar conclusions have been made in [22] for a model of non- linear massive gravity in Stuckelberg formalism. 
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5 Conclusions 



In this paper we have presented a chiral model for ghost-free massive (bi-)gravity. The theory we 
considered is given by two copies of the Plebanski action for general relativity, plus an interaction 
term. On shell, the interaction corresponds to the sum of two terms (|22|) . one of which coincides with 
the symmetric interaction between tetrads considered in [25] (second item of ([2D) and the other features 
an alternative gauge-invariant contraction of the tetrads. We performed a complete canonical analysis 
and proved that the theory only propagates seven degrees of freedom, hence effectively removing 
the Boulware-Deser ghost usually plaguing models of massive gravity. Thus, our results give an 
independent and explicit proof of the absence of the ghost for a special type of interaction terms 
O [T2| I13j . and furthermore extend it for the first time to a first order action of gravity, with 
independent tetrad and connection variables. 

It is instructive to compare and contrast the formulation of bi-metric gravity used here with the 
other, more standard ones such as metric and tetrad formulations. The former is the most economic one 
in terms of fields, as any gauge formulation introduces unphysical redundancies, and both the tetrad 
and chiral Plebanski-type formulations have more of these. However, the interaction terms needed to 
remove the ghost are non-polynomial in the metric framework, which makes computations increasingly 
complex. The tetrad formulation converts these complicated interactions into simple polynomials, but 
at the same time it introduces the additional gauge symmetry of local frame rotations without really 
simplifying the Hamiltonian analysis. What we have considered in this paper is a chiral model, which 
lies in between the pure metric and the tetrad formulations since it brings only 3 local gauge rotations 
as compared to 6 in the tetrad formalism. The reduction of the gauge group tremendously simplifies 
the canonical analysis, a fact that can be exploited to perform a number of explicit calculations. 

Our formalism does not allow us to incorporate the first and third terms in ([2]), hence it may look 
more restricting, allowing de facto only the second term in the list. However, as we discussed, the 
use of 2-forms as fundamental variables also opens the door to the possible existence of an infinite 
class of ghost-free theories, by turning the cosmological constants into arbitrary functionals of the 
Lagrange multipliers ^ ,=t . These in turn introduce new interactions for the metrics, and the procedure 
can be argued to not introduce extra degrees of freedom. The existence of a more general class of 
ghost-free interaction is interesting, and deserves to be further studied. 

The chiral formulation has a drawback in that the analysis of bi-metric theory holds, strictly 
speaking, only for metrics with Euclidean signature. The calculations extend immediately to the 
physical case of Lorentzian signature, however in this case one works from the start with complex 
fields, to be supplemented by appropriate reality conditions. The simplest possibility is to borrow the 
usual reality conditions for Lorentzian Plebanski gravity and apply them independently to both sets 
of fields. This almost works: it induces real and Lorentzian metrics, and the reality conditions are 
preserved by time evolution. However, the total action is not made real by these conditions, unlike in 
the simple gravity case. One imaginary term survives and spoils the theory, and this is precisely the 
parity-odd extra term in (I22D . Interestingly, the term that makes the standard reality conditions fail 
is the same that characterizes the difference between bi-gravity in the metric and tetrad formulations. 
Hence, finding more non-trivial reality conditions could also be related to finding a reformulation of 
the Plebanski version totally equivalent to the metric one. 

A way to avoid the problem of reality conditions would be to start with the non-chiral Plebanski 
theory based on the full Lorentz group [60\ [4"3"1 [4"4"] , which is at the basis of the spin foam quantization 
program to general relativity [451 146] . However, its canonical analysis is of the same difficulty as the 
one in the tetrad formulation, see \47\ I48j . Furthermore, modifying the non-chiral Plebanski theory 
in a way similar to the one considered in section 14.61 in contrast to the chiral case, one obtains a 
theory with generically 8 degrees of freedom [49], that can be interpreted as those of a generic massive 
gravity including the scalar ghost [50] . This suggests an interesting program of realizing the ghost-free 
massive gravity as a particular subclass of modified (non-chiral) Plebanski theories. An investigation 
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along these lines has appeared in [51] . 

Another interesting open question concerns the existence of stable subsectors with fewer degrees 
of freedom. Apart from the trivial cases of a = 0, oo (corresponding respectively to two decoupled 
gravity theories and a single one), our analysis allows us to explicitly characterize such subsectors by 
analyzing the stability condition for the secondary constraint. They appear unstable at first sight, but 
a more complete study of their dynamical properties and geometric interpretation is needed to arrive 
at a conclusive result. 

Finally, it would be interesting to see whether quantization techniques, such as loop quantization 
and spin foam approach, developed to deal with gauge formulations of gravity, can be extended to the 
models describing massive (bi-)gravity. 
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Appendix 

A Evaluation of the secondary constraint 

In this appendix we evaluate the secondary constraint whose general form is given in (|4ip . The first 
term in this expression is just the Poisson bracket between the diagonal and off-diagonal Hamiltonian 
constraints computed in ()36c|) . It can be further simplified as 



{%,%} = -4a 



e[Ef- H°E J b H° )+h[Hf- E 4 Hi El 



pa rrj jpb 



3 rib ^1 



(53) 



where e = det-E", h = detHf 1 , and thus reproduces the first term in (|42p . The second term includes 
contributions of the remaining three terms in (1411) . To evaluate them, the first step is to invert the 
matrix 

eEl - hHl 



A 1 

a 



±{A,Gi} = ^e flbc (E b 3 Ei 



Tjb TTC 

H j H k 



Its inverse can be computed as follows 
where 



1 

2A 



eEf + hH° 
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fpa fpb jpb jpa\ r jkl r ttc ird 

Ei Ej - &{E i j e e bcd H k H l 



A = det A l a = (e - hT)e - (h - eT)h, 
and T = H l a E?, T = E l a Hf. Using this result, one further computes 



A' 1 
A' 1 



e (E*E$ + H?H% )+{h- eT)H?Et e^E^Hf, 



ra Tjb 



a rpb 



HC Tjd 



h[E?E\ + H*Ht)+{e-hT)H?E 



TO, Tjb 



ra rpb 



C TTCl 

3 3 ' 



(54) 

(55) 
(56) 

(57) 



which allows to get the last term in (|41|) in the following form 

{H, G 3 }{G 3 ,V b }- l {V b , V a }{A, GiY^Guii} 

= 4aA"V (Ejd s E! - Hld s H?) e, cd + H«h() E c k H d k ej pq H?E( ' EP m H^. 



(58) 
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This contribution can be combined with ^-dependent terms from 



-4aA" 



e {EtE\ + HfHn +(h- ef)HtE\) A^H 



h (ElEl + H?H\) + (e - hT)H?E*) A^jj e^E^H^fEf H[ 
- terms linear in 77. 



(59) 



One can check that these terms complete the derivatives in (|58p to covariant ones, which in turn can 
be replaced by the two Gauss constraints provided one subtracts the missing terms linear in r\. The 
contributions proportional to the Gauss constraints can be dropped, whereas the terms linear in r\ are 
combined with similar ones from (|59|) producing the last term in (|42|) explicitly given in 



B Perturbative analysis around H a = 

In this appendix we analyze in detail the secondary constraint and its stability condition in a first few 
orders of a perturbative expansion around H a = 0. It should be emphasized that this analysis is more 
general than a linearization around a flat background. In particular, we will show that, in contrast 
to the later (under the assumption of non-degeneracy of the physical metric), it allows for solutions 
of the conditions (|52p for which the stability condition of the secondary constraint does not seem to 
impose a restriction on the Lagrange multipliers. 

First, let us compute the secondary constraint up to quadratic terms in H a . Taking into account 
that 

A = e 2 + 0(H 2 ), (60) 
from the explicit expression ()42|) together with (|44p one can show that the constraint can be written 



as 



* « — 4ae 



;i + 0(H 2 )) E*rf a + vlHfEiH] - rf a H*& b Hf - r,lH?E?H>(EE) bc + 0(H 3 ), (61) 



where the constraint Qi has been extensively used. As a result, dividing by the prefactor of the first 
term, one concludes that in our approximation it is sufficient to consider 

* = Efrtl + * (2) , * (2) = VIH?EIH* - rfM^R* - rj^HfEjHj(EE) bc . (62) 

In particular, we observe that at the leading order the secondary constraint ensures the vanishing of 
the trace part of the r\ field. Similarly, using (134^ . one can interpret the "off-diagonal" Hamiltonian 
constraint % as a restriction on the trace part of H (since the last term gives aE^Hf), whereas Qi and 
T> a can be viewed as restrictions on the antisymmetric parts of r] and H, respectively. This observation 
is already sufficient to conclude that in the (H, rj) sector we can have at most 5 propagating degrees 
of freedom, contained in the symmetric traceless parts of the fields, which are expected to correspond 
to the degrees of freedom of a massive graviton. The ghost scalar carried by the trace part is removed 
by the constraints % and ^> . 

Next, let us analyze the stability condition of the secondary constraint which requires the vanishing 
of the Poisson bracket {^f,H to t}- We will compute it keeping only terms which are at most linear in 
H a . To this end, using ([ID]) and (|57p. on the constraint surface we obtain 



G a = - e^NEfEfep^Hf + 0(H 2 ), 

j,i*N (D a H? - e-\ abc A^EfE*E b k m) + e^Ge^EfE^ {h^ - HfEg) + 0(H 2 ). 
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Given the above results, a direct computation gives 



+ 



8(H N + H G ) dV® 



- + 0(H' 



N 



+ 



2e tjk EfH*e klm V l a r ] F + Sae^^E^E^H^ + e^EfE^ 1 " 1 ^ 



brf a dH: t 



a ipb ttc i ijk fpa rpb Mm m 



+0 a Ge^it? m^Hfri - H V k 



+ 0(H 2 ). 



(64) 



Using (|62p and the constraint % in the second term, on the constraint surface (including this can 
be further simplified as 



(65) 



fftot} « # ((6ae - ^r?*^V a ) £ c fc # fe c + 2Ef4H^ a 
+ G (6ae - + W H ) + 0(# 2 ), 



where 



e« fc ^? U d H?E k D a E b - E'MHfDa 



Vb 



+& d HfriD a El + 2E b J EU d HfE™D a E^ + Ep^^Hf) + D a ( V k E?H?Ei 



(66) 



In particular, all terms with derivatives of the Lagrange multipliers cancel in agreement with the 
general conclusion of section 14.41 

Both coefficients, in front of N and G, appear to be non-trivial functions on the phase space. If 
at least one of them is non- vanishing, the equation {fy,H tot } = can be considered as an equation 
for the corresponding Lagrange multiplier. Then the Dirac's procedure is completed and ^ together 
with one of the Hamiltonian constraints is of second class. 

Let us now analyze the situation where both coefficients in (|65p simultaneously vanish. One might 
expect that in this case the constraint structure changes drastically. In particular, there is a chance 
that both Hamiltonian constraints and become first class. In such situation, although the ghost 
degree of freedom is still removed, it would happen due to an additional gauge symmetry, rather 
than by means of constraints which are of second class. The perturbative expansion presented in this 
appendix provides a framework to explicitly address these issues. 

First, we need to understand whether it is really possible to have both coefficients vanishing. 
Explicitly this requires 



y := ( 6ae - E?4&ivL) E k H c k + lE^rfM^ = 0, 



y := 6ae - E?4E b r? a + W H = 0. 



(67) 



To start with, we consider a linearization around a flat background which can be described by Ef = 8f 
and A a = 0. Using the flatness of the connection, from the Hamiltonian constraint % (|34p one finds 
that 

Etrf h E)rf a = 0(H). (68) 

Combined with y in (j67[) . this leads to 6a = 0(H) which is inconsistent with our assumption of H a 
being infinitesimally small. Thus, near a flat background y ^ so that % forms a second class system 
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with the secondary constraint and the model describes one massless and one massive graviton, 
consistently with previous findings [25j . 

The above consideration suggests that one should look for a common solution to (|67|) and all other 
constraints in the sector with a constant curvature set by the mass parameter a, rather than the flat 
one. Let us do this setting for simplicity H a = which automatically ensures y = 0. In the sector 
with vanishing H a the constraints become 

¥ = Efni, Qi = D a E?, Qi = e ijk r&E ak , 

C a = E*>(F l ab + ^ k rf a r, k b ), C a = E*(D aV l - D b rf a ), (69) 
H = g e E i E j( F ab + e VaVb ), H = e J E t EjD a ri b . 

Introducing rfi = rf a E c -, and requiring that it is symmetric and traceless, solves Qi and \& so that we 
remain with the following equations 

l ahc ^ a Fl = 6a, D a E? = 0, 

E\Fi b = 0, Etfi (D a g - D b E{) - E\W a D b ^ = 0, (70) 



rf 3 r]ij = Da, rj 

The first two of these equations require the connection to be consistent with the triad and to have a 
constant curvature, whereas the fifth equation fixes the "norm" of the field rf 3 . Then we remain with 
5 differential equations on 9 components of Ef. It is natural to expect that they can be simultaneously 
solved which shows that the space of common solutions to the conditions (|67p (or more generally to 
(I52p ) and the constraints is not empty. 

The next step is to study the stability of the additional constraints ()67|) . Our approximation allows 
to compute the commutators in the zeroth order in H a . Explicitly, one finds 



6a, r] ij ~e abc EiD a El = 0. 



{y, fltat} ~ {Hn + n G ,y} + 0{H) « -Wrftfrfi + 2Gii k E*E) m 4 m D a E£ + O(H), 
{y, fltot} ~ {Hn + Ha, y} + 0(H) « 2Ne i i k E?E b jVkl r ] lm D a E$ l (71) 



+ G 



2ii k EU 3 E\F k b - || A z y^ h D a {Et^)^W H 



+ 0(H), 



All terms here appear to be non-vanishing. As a result, the stability of y and y cannot be achieved 
unless both Lagrange multipliers N and G simultaneously vanish. However, we are not interested in 
such situation as it corresponds to degenerate metrics. Therefore, the sector defined by (j67]) does not 
seem to possess an interesting dynamics and probably does not have any physical significance!^! 



C Commutator of two Hamiltonian constraints 

Let us consider the commutator of the two first class constraints playing the role of the Hamiltonian 
constraint. Such constraint can be obtained by collecting all terms in the total Hamiltonian (|38p which 

10 An alternative possibility to treat the conditions (|67[) would be to insert them into the original Lagrangian. Then 
they appear in the total Hamiltonian with appropriate Lagrange multipliers and their stability should be studied at the 
same step as the stability of the primary constraints. In this scenario, the stability of y and y is ensured by fixing the 
Lagrange multipliers and does not require the vanishing of TV and G. However, due to 

{y(x), y(y)} w W j rf k r 1 k i P(x - y) + 0(H), 

the constraints y and y appear to be second class and affect the expression of the secondary constraint. As a result, this 
scenario leads to a constraint surface which does not intersect with the one corresponding to the physically interesting 
case of massive gravity. 
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are proportional to the "diagonal" lapse function N after plugging in all solutions for the Lagrange 
multipliers fixed by the stability conditions. Thus, it is given by 



-j/f.C. 



(72) 



where G(N), G a (N), ip l (N) are those solutions which can be found from (|4L)p and (|51|) . Our aim here 
is to get the function analogous to N a (N,M) in (|19p to read off the physical metric determined by 
the diffeomorphism algebra. Therefore, we are not interested in the full commutator, but only in the 
terms proportional to the diffeomorphism constraint T>. Using the constraint algebra presented in 
(f36l) . one finds 

= V a (N u N 2 )V a + • • • , (73) 

where 



V B (J\Ti,JV 2 ) =V a (N u N 2 ) + F a (G(iV 1 ),G(^ 2 )) - U a (N 1} G(N 2 )) - U a (G(JVi),J\T 2 ) 
+ L a (G(N 1 ),G(N 2 )). 

Let us take 

G a (N) = Nf a + G(N)g a , G(N)=KN, 

where our results imply that 

f a = {n,0i}{&,i> a }~ 1 , g a = {n,Gi}{Gi,f>a}-\ n = -y/y. 

Then one easily calculates that 

V a (N u N 2 ) = K ab (N x d h N 2 - N 2 d h N x ) , 

where 



K ab = (1 + TV) [{EE) ab + (HH) ab j - 2TZ lE?H° + E b H?) + (f a + lZg a ){f + Tlg b ) 
= (Ef - KH?)(E* - UH b ) + {Hf - nEf)(H b - UE b ) + (f a + Kg a )(f b + Kg b ). 



(74) 

(75) 
(76) 

(77) 
(78) 



The function K, ab is expected to encode the spatial part of the physical metric of the coupled theory 
as K, ab = gg ab . The result (|78ft clearly shows that the physical metric is a very complicated function 
being expressed in terms of "diagonal" and "off-diagonal" variables. 

Using results from appendix [Bj it is possible to get a more explicit representation for K, ab in the 
quadratic approximation in H a . One finds 



ab 



(EE) 



ab 



1+ + 



2E9r^Hpj 



c 



6ae- E 9 k tfE[r] l g/ 



(HH) 



ab 



o | rpo, Tjb _j_ rpb rra 



2Efr ] i d Hfrf c 



(79) 



v 6ae - E 9 k r] k f E(r] l g 

+ e ijk E^H^E^e lmn E^HiE^ + 0{H 3 ). 

Expanding around a bi-flat background 

% = % + f~A~H~ v ~o(l), 

this further simplifies to 

K ab ^ § ab + f a + f a + f a f b + fiftyb 

+ 5 ab {Htf + 2{H% + H b a )Ht + e ajk H k e bmn H^. 



(80) 



(81) 
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The expression in the first line gives the fluctuation of the metric defined as a symmetric combination 
of the two metrics constructed from and expanded around a bi-flat background. All terms in 
the second line contain contributions either from trace of H or from its antisymmetric part. These 
components are constrained by the constraints % and T>, respectively. However, they are restricted 
not to vanish, but to be given (in our approximation) by derivatives of n: 

Hj = ± ^ k d^, e ajk H$ = --L di rf a - (82) 

Thus, we found that although in the linear approximation the fluctuations of the diagonal triad Ef 
are consistent with the fluctuations of the physical metric, already at the quadratic order the metric 
gets contributions from all fields. 
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